Logrithmic corrections to the RG flow for the two-dimensional bond disordered Ising 

model 
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Using the mapping of the partition function of the two-dimensional Ising model onto a pfaffian we 
evaluate the domain wall free energy difference for the pure and disordered Ising model close to the 
pure fixed point. Using this method very large lattices can be studied exactly and we confirm that 
disorder even including frustrating interactions indeed are irrelevant close to the pure fixed point. 
The finite-size renormalization group flow shows a power-law behavior modified by a logarithmic 
term that dominates for small lattice sizes. 
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I. INTRODUCTION 

The influence of disorder on the pure two-dimensional 
Ising model has been a problem of long standing interest. 
The works of Ludwig [1] , Shankar [2] and Dotsenko and 
Dotsenko [3] have shown that logarithmic corrections to 
the specific heat arises once disorder is present but the 
otherwise the critical exponents remain unchanged from 
the pure case. In a certain sense one might say that dis- 
order is marginally irrelevant for the pure Ising model. 
From the Harris criterion [4] we know that disorder is rel- 
evant if the correlation length exponent, v > 2/d. How- 
ever, for the two-dimensional Ising model v — 1 and the 
Harris criterion gives no definite answer to relevance or 
irrelevance of disorder which in this case is a marginal 
perturbation. Most of the analytical treatments of this 
problem considers only non-frustrating disorder and one 
might ask the question if frustration could have a differ- 
ent influence on the system than non-frustrating disorder. 

In this paper we revisit this problem and by mapping 
the partition function onto a pfaffian we are able to ob- 
tain large scale numerically exact results independent of 
the disorder distribution for various different quantities 
close the pure fixed point. To our knowledge pfaffian 
techniques have so far not been applied to fully disor- 
dered systems. Compared to more standard transfer ma- 
tric techniques this method presents significant advan- 
tages notably in the system sizes it is possible to treat. 
In some cases it is possible to obtain exact results for 
system sizes comparable to what has been treated with 
Monte Carlo methods. Using finite-size scaling ideas we 
show that disorder indeed is irrelevant at the pure fixed 
point even including frustration. However, the finite-size 
renormalization group flow has a novel logarithmic be- 
havior recently predicted to be rather universal for many 
models by Cardy [5] . For small system sizes disorder ap- 
pears to be relevant and only for very large system sizes 
does it become apparent that the flow is towards the pure 
fixed point. 



II. METHOD 

Our approach is inspired by the domain wall renormal- 
ization group (DWRG) developed by McMillan [6] and a 
recent modification of this, used to determine the fixed 
point structure of the two-dimensional disordered Potts 
model [7]. In this approach the stiffness of the system is 
calculated and a finite size scaling ansatz is used to deter- 
mine the flow. We briefly recapitulate the central results. 
In the absence of disorder it is known from hyperscaling 
that the singular part of the free energy density scales as 
the inverse of a correlation volume 
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Exploting the point of view of the DWRG we assume 
that the critical part of the free energy, f s , comes from 
the difference in free energy between two configura- 
tions with different boundary conditions f s = Af = 
/periodic - /antiperiodic- Using standard finite-size scaling 
arguments [7] we find that the total free energy difference 
AF scales as 
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the 



with g a universal function and <5 = \T — T c \. At T c 
free enrgy difference is thus a universal amplitude. 

Ag{8L 1 / u ) can be regarded as the stiffness of the sys- 
tem, p{L,T). In an ordered phase we would expext p 
to diverge with the system size, whereas in a disordered 
phase p should scale to zero with the system size. Deep in 
the ordered phase it is known that the leading correction 
to AF diverges linearly with the system size L [8]. 



\AF/T\ ~ L 



In 



1 1-1*1 
\z\l + \z\ 



tanh ( — 
T 



(3) 



At T — T c this term is rigourously zero and the finite 
universal free energy difference comes from higher order 
corrections. The critical point, T c , can now be located by 
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standard methods, i.e. by tracing AF as a function of T 
for several different system sizes and locating the point 
where the lines cross. 

The Hamiltonian that we consider is the bond disor- 
dered two dimensional Ising model: 

H = - J2 J V S i S i> Si=±l. (4) 

<i,3> 

To facilitate some of our calculations we restrict our 
treatment to a bimodal distribution of the bonds: 

V{J)=x8{J -l) + {l-x)S(J -1), (5) 

where x is the concentration of antiferromagnetic bonds. 
However, one should note that the present techniques 
can be used for any distribution of disorder, continuous 
as well as discrete. 
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figure 1: The antiferromagnetic seam on the torus. 

Since the works of Kac and Ward [9] it has been known 
that the partition function for some two-dimensional 
models can be expressed as a pfaffian. A general the- 
orem due to Kasteleyn [10] shows that this is possible 
for all planar graphs. For the problem at interest here 
this means that a pfaffian formulation is not possible in 
higher dimensions. The Pfaffian formulation for the two- 
dimensional Ising model was later developed in detail by 
several authors [9,11,12,8]. We refer the reader to the 
litterature for the derivation of this mapping. What is 
very noteworthy is that the the mapping onto a pfaffian 
can be exploited without problems even for the fully dis- 
ordered case. For the disordered Ising model on a torus 
one finds a sum of four pfaffians: 

Z = 2 i2 (cosh(J/T)) 2i2 - Ar (cosh(-J/T)) JV 

(-Pfi + Ph + Pfa + PU)/2, (6) 

with N the number of AF bonds. In order for this expres- 
sion to be useful it is necessary to be able to evaluate the 
pfaffians for rather large matrices. Usually this is done 
by using the well known relation for a skew-symmetric 



matrix D, \Pf(D)\ = Det(D) 1 / 2 . However, even though 
efficient methods exists for evaluating determinants, this 
expression only determines the absolute value of the pfaf- 
fian. Fortunately, recent developments in combinatorial 
algorithms have shown that it is possible to evaluate the 
pfaffian directly in polynomial time [13], without mak- 
ing use of the above formula. This polynomial algorithm 
is relatively slow compared to the more efficient methods 
for evaluating the determinant and for the larger systems 
sizes we have considered we have calculated the determi- 
nant and obtained the sign using continuity arguments. 
This can be done since the sign of the four pfaffians is 
positive at sufficiently high temperatures. 

III. RESULTS - PURE SYSTEM 

As an illustration of how well this works we first con- 
sider the pure Ising model without disorder and try to lo- 
cate the known critical point, T c = 2/ log(l + V%), using 
the above mentioned techniques. Introducing a line of an- 
tiferromagnetic bonds along one line through the system, 
as illustrated in Fig. 1, we can calculate the total free en- 
ergy difference, AF, between periodic and antiperiodic 
boundary conditions. Our results are shown in Fig. 2 for 




figure 2: |AF/T| for L = 16,24,32,40,48 as a function 
of T. 

system sizes of L = 16,24,32,40,48 as a function of T. 
One should note that the scale on the x-axis which is ex- 
tremely close to the critical point, T c = 2.269185 .... On 
this very fine scale it is clear that corrections to scaling 
are present. Only in the thermodynamic limit do we ap- 
proach the universal value of the stiffness at the critical 
point. However, even neglecting these corrections, the 
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critical point can be located to four decimal plase. The 
slope of AF at T c should scale as L 1 / 1 ' and can be used 
for a very precise determination of the exponent v. Since 
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figure 3: |AF/T| at T = T c for L = 16, 24, 32, 40, 48. The 
open circles (o) indicate the exact numerical results and 
the solid circles (•) indicate a power-law fit to the data 
of the form 0.986485 - 0.303151£- 2 ' 00853 . 

T c in this case is known we can explicitly study the cor- 
rections to scaling at T c by simply calculating AFt=t c - 
Our results are shown in Fig. 3. A simple power law fit 
to the exact numerical results works extremely well and 
we find that AFt=t c approaches a universal constant in 
the following manner: 

\AF/T\ ~ 0.986485 - 0.303151L" 2 00853 . (7) 

From conformal invariance it is known that for the geom- 
etry of an infinite long cylinder AF/T = wr] [14], where 
r\ is the magnetic critical exponent. Hence, for this ge- 
ometry, AF directly measures a bulk critical exponent. 
However, for the toroidal geometry used here we are not 
aware of any results for the universal number, 0.986485, 
even though many exact properties are known for the 
Ising model. We would welcome any information as to 
where and if such a calculation has been performed. The 
corrections to scaling follow a standard power-law form 
with an exponent extremely close to 2 as expected. 



IV. RESULTS - DISORDERED SYSTEM 

We now turn to a discussion of our results in the pres- 
ence of disorder. In this case we need to consider not 
only AF but also its standard deviation <r(AF). The 
relevant variable to consider is then r = a(AF)/[AF] 
where [•] denotes the disorder average. We are interested 
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figure 4: The slope at T = T c for L = 16,24,32,40,48. 
The open circles (o) indicate the exact numerical results 
and the solid circles (•) indicate a power-law fit to the 
data of the form 2.941og(L) - 36 L" 070 . In the inset we 
show the same results as well as the extrapolated form 
out to L = 10 3 . 

in studying the behavior of this quantity close to the pure 
fixed point where a(AF) is trivially zero. In order to de- 
termine the relevance or irrelevance of the disorder we 
need to calculate: 



slope = 



d (j(AF)i 
dx [AF] 



£C=0- 
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Here x denotes the concentration of antiferromagnetic 
bonds. At a standard fixed point we would expect the 
slope to follow a power-law behavior L A and in the case 
of A < we qualify the perturbation as irrelevant and as 
relevant if A > 0. The calculation of the derivative with 
respect to x is performed by introducing a single anti- 
ferromagnetic bond into the lattice. For the case of the 
bimodal distribution that we consider it turns out that 
there are only two non-equivalent positions on the torus. 
Hence, we can evaluate both [AF] and a(AF) exactly. 
A single antiferromagentic bond corresponds to a given 
value of x and since r at x — is trivially zero we can 
numerically determine the derivative. One should note 
that the error introduced by numerically determining the 
derivative is extremely small due to the large lattice sizes 
and the corresponding tiny value of x. 

Our results for the slope at T = T c are shown in Fig. 4 
for lattice sizes of L — 16, 24, 32, 40, 48. Surprisingly the 
slope increases with the lattice size L and one could be led 
to the conclusion that in the presence of this frustrating 
disorder the perturbation is marginally relevant instead 
of marginally irrelevant as is known to be the case for 
unfrustrating disorder. However, since our results are 
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numerically exact we can verify that the functional form 
is not a simple power-law. Instead, we find that the only 
acceptable form is: 

AlogiLYL-^ (9) 

with a ~ 0.36, (3 ~ -0.07 and A = 2.94. This is 
a functional form quite well-known for the correlation 
functions [1,2] in such disordered systems and as re- 
cently shown by Cardy [5] it is common for many low- 
dimensional models. Considering the precision of our nu- 
merical results we believe this form to be correct up to 
small errors in the exponents. If this is indeed the case 
the 'flow' will not reverse until lattice sizes of approxi- 
mately L=200 are reached, a system size not accessible 
by our methods. We do not attach any importance to 
this length scale since it depends on the form of the dis- 
order and various other microscopic parameters. In the 
inset in Fig. 4 we show our numerical results along with 
the extrapolation out to lattice sizes of L= 10000. 
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V. CONCLUSION 



Applying the well-know mapping of the partition func- 
tion of the Ising model onto a pfaffian we have succeeded 
in evaluating the finite-size renormalization group flow 
around the pure fixed point in an almost exact numeri- 
cal manner. We have used a bimodal distribution of the 
disorder which includes frustration for which the domain 
wall free energy as well as its standard deviation can be 
evaluated exactly. Our results show that very large lat- 
tice sizes have to be used before it becomes clear that 
disorder, even frustrating, is irrelevant. It is possible 
that the functional form for the finite-size renormaliza- 
tion group flow Eq. (9) can be verified using analytical 
methods. 
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